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Based on exactly mapping of a many-body electron-phonon interaction problem onto a one-body 
problem, we apply the well-established nonequilibrium Green function technique to solve the time- 
dependent phonon-assisted tunneling at low temperature through a single-molecular quantum dot 
connected to two leads, which is subject to a microwave irradiation field. It is found that in the 
presence of the electron-phonon interaction and the microwave irradiation field, the time-average 
transmission and the nonlinear differential conductance display additional peaks due to pure photon 
absorption or emission processes and photon-absorption-assisted phonon emission processes. The 
variation of the time-average current with frequency of the microwave irradiation field is also studied. 
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I. INTRODUCTION 

The role of inelastic scattering in electron transport 
through mesoscopic devices is an active topic of current 
theoretical and experimental interest .LSiiMiSiSiLSiiliilLii 
Thanks to the modern advanced nanostructure tech- 
niques, we can build semiconductor quantum dots (QDs) 
with easy control of its size and shape. Recent experi- 
ments found that the electron-phonon interaction (EPI) 
for longitudinal optical phonons in semiconductor QD 
can be substantially enhanced due to the quantum con- 
finement effect, and also revealed the importance of the 
EPI in transport measurements through this kind of sin- 
gle electron transistor. On the other hand, a growing 
interest has been payed recently to electrical transport 
through a very small single molecular device. In con- 
trast to the semiconductor QD, the molecular material 
possess much weaker elastic parameters, leading to that 
their internal vibration degrees of freedom is of "low- 
energy", 1-10 meV i 1 ' 2 ' 3 ' - Therefore, when there are elec- 
trons incident into the molecule through tunneling junc- 
tion from external environment, these low-energy boson 
modes (phonons) can be easily excited due to the strong 
coupling between these modes and the molecular elec- 
tronic states, and thus react inevitably back to the tun- 
neling of electrons even at very low temperature. For 
example, the interesting experimental results reported in 
Ref.Q demonstrated that electron transport through a 
single Ceo molecule was significantly influenced by a sin- 
gle vibrational mode. 

This phenomena has provoked a large amount of theo- 
retical work on the problem of tunneling through a single 
level with coupling to phonon modes based on the Fermi 
golden rule^ the kinetic equation approach*^ and the 
nonequilibrium Green's function techniquesAiSiii Natu- 
rally, both the phonon energy and the electron-phonon 
coupling strength are of importance to determine trans- 
port properties of these single electron transistors. Most 
theoretical attempts to study this transport problem so 
far have focused on stationary case. As far as the energy 



of boson in molecule is considered to be in the region 
of GHz to THz, one can imagine that more rich physics 
could be exploited if the device is subject to a microwave 
(MW) irradiation field. It has been reported that the 
microwave spectroscopy is a possible tool to probe the 
energy spectrum of small quantum systems and to mea- 
sure the decoherence time of the quantum states^ So 
the photon-assisted tunneling (PAT) could provide a new 
way of understanding more essence of the EPI influence 
on the transport properties in molecule. 

The purpose of this paper is to study the time- 
dependent coherent transport through a single-molecular 
device coupled to a local dispersionless phonon mode by 
means of nonequilibrium Green's function (NGF) under 
the adiabatic approximation^ Acutally, when inelas- 
tic processes such as phonon emission and absorption 
are considered in the tunneling event, electron scatter- 
ing becomes a many-body problem involving electrons 
and various excited phonon states of the system. Since 
the coupling between electrons and phonons is strong 
in the semiconductor and the single-molccuar QD, the 
usual perturbation theory is invalid for dealing with this 
scattering problem, although it has been extensively ex- 
ploited with enormous success to study the conventional 
transport in bulk material. Recently, a nonperturba- 
tive scheme has been proposed by Bonca and Trugman 14 
for studying the small polaron of the coupled electron- 
phonon system described by the Holstein modeU^ and 
the Su-Schrieffer-Heeger (SSH) models The fundamen- 
tal idea of this method is to rewrite the Hamiltonian in 
terms of the combined electron-phonon Fock space, such 
can exactly map the many-body problem onto an one- 
body scattering problem. Later on, this method has been 
further applied to the topic of inelastic electron scatter- 
ing in mesoscopic quantum transport in semiconductor 
QDs and molecular wires under the framework of the 
Landauer-Buttiker scattering theoryi 1 i 1 £ In this paper, 
we will redirect this procedure in terms of the NGF tech- 
nique and address the inelastic PAT in a molecular QD. 

The rest of this work is organized as follows. In section 
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II, we introduce the Hamiltonian for resonant tunneling, 
where electrons interact with the boson fields localized on 
the QD or molecule. Then we describe our methodology 
for studying the time-dependent tunneling and comput- 
ing the time-average transmission and current by means 
of NGF technique. In section III, we give the numerical 
results and discussions. We find some new features in 
the time-average transmission spectrum associated with 
the combination effects of the photon-phonon emission or 
absorption. Then we predict that this photon-phonon- 
assisted tunneling can be observe experimentally in the 
bias voltage-dependent differential conductance and in 
the microwave spectroscopy of dc current. Finally, a brief 
summary is presented in Sec. IV. 



II. METHOD AND FORMULATION 

A. Model and Hamiltonian 

We consider the system under investigation as the sim- 
plest case: a single-site QD couples to two noninteracting 
reservoirs via tunneling and interacts with a dispersion- 
less optical phonon localized in this site. It is anticipated 
that this is sufficient to illustrate the main physics for 
transport through the semiconductor QD and molecule 
in the presence of many boson modes. Moreover, we 
neglect the spin degree of freedom and any effects of 
electron-electron Coulomb interactions. Therefore, the 
Hamiltonian of this system can be split into three parts, 
(1) the two isolated leads H\ ea a, (2) the single-site QD 
Hqd, and (3) the tunneling part Ht\: 

Hqb = t d {t)c\c d + hw ph tfb -\c\c d {b +fet),(lb) 

= E^(4 Cd + H - c -)- ( lc ) 

r),k 

c \)k (°vk)i c d (cd) are the creation (annihilation) operators 
for the noninteracting electrons with momentum k in the 
r] (= L/R) lead, and for the electronic state at the single- 
site, respectively. The role of the MW fields irradiated 
on the whole system can be described by a rigid shift of 
the single-electron energy spectrum under the adiabatic 
approximation^ e v k(t) = e" fe + v v (t) (77 = L, R) and 
ed(t) — €d + Vd{i), e° fe and td are the time-independent 
single electron energy without MW fields, v n /d{t) is the 
time-dependent part, v v /d(t) = v v /d cos fit, with v v /d be- 
ing the irradiation strength in different elements of the 
device and f2 being the frequency of the MW field. Huj p h 
is the energy of the dispersionless phonon mode and A is 
the on-site EPI constant. In the present work, we neglect 
the direct coupling between the MW field and the boson 
field. V v stands for the tunneling coupling between the 
dot and 77th lead. 



Obviously, this problem described by the above Hamil- 
tonian Q is a many-body problem involving the phonon 
emission and absorption when the electron tunnels 
through the central region. Following the pioneer work 
of Bonca and Trugman, we can expand the electron 
states in the single-site onto the polaron eigenstates, the 
direct-product states of the single-electron states and the 
phonon Fock states, 

|0,n) = ctffio>, (n>0) (2) 
v nl 

which means that the electron on the site is accompa- 
nied by n number of the phonon hu p h (|0) is the vacuum 
state). After performing this representation, the many- 
body on-site Hamiltonian (|Ib(l can be mapped onto an 
one-body oneti^iiiS 

Hqb = E [fa + nhu p h)\0, n) (0, n\ 

n>0 

-AVn+T(|0, n + 1) (0, n\ + |0, n)(0, n + 1|)]. (3) 

Similar expansion of electron states in the two non- 
interacting leads with the combined space \rjk,n) = 

c^ k ^-L. |0) will change the simple 77 lead Hamiltonian 
Qlafl to a pseudo-multi-channel model labeled by the 
phonon quanta n with a weight factor P n = (1 — 
e -t3w ph ^ e -nf3u ph ^ which is the statistic probability of the 
phonon number state \n) at the finite temperature T 
(}3 = 1/k^T). Notice that the formula ^ n P n — 1 guar- 
antees the statistic properties of the leads unchanged (see 
the details in the following subsection). 

Finally, the tunneling part Eq. (fTcf> can be also rewrit- 
ten in terms of this basis set. Considering the presump- 
tion presented here that there is no exchange between 
electrons and phonons in the tunneling process, only hop- 
ping of the electron states are involved in the effective 
tunneling Hamiltonian: 

Hti = Yl V r rjn (|»?fe s »)(0,n| +H.c), (4) 

77, A:, n 

where l^fc) denotes the electron state with momentum k 
in the 77 lead. It should be noted that this presumption 
is rational because the high-order tunneling processes ac- 
companied by the phonon emission and absorption are 
much weaker than the direct tunneling events. And in 
fact, inclusion of these high-order processes in the calcu- 
lations is an instant task under the present theoretical 
framework. V nn is the coupling between the nth pseudo- 
channel in the 77 lead and the QD. 

We depict the many-body tunneling process after re- 
modeling EPI system in terms of the combined electron- 
phonon Fock space: An electron incident in the pseudo- 
channel labeled by a particular phonon number n in the 
left lead can transport elastically to the central region, 
then excites or absorbs 771 numbers of phonons and ex- 
its inelastically into the (71 ± m)th channels in the left 
and right leads; or experiences no exchange with phonon 
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FIG. 1: Schematic description of the inelastic scattering prob- 
lem for a single site with on-site EPI. Each phonon state of 
the site along with the Bloch state of the electron in the lead 
can be visualized as a pseudo-channel labeled by n, which 
is connected to the two leads with the hopping parameters 
Vl/r m the wide band limit. These channels are connected 
vertically by the EPI A. 

and leaves elastically at the same channel in both leads. 
Fig. 1 shows a graphical illustration for this description. 

So far, all theoretical solutions in the literature for the 
phonon- assisted tunneling under above mapping Hamil- 
tonian are based on the Landauer-Biittiker scattering 
formulation^iiii& In this paper, in order to utilize 
the well-established NGF technique to deal with the 
time-dependent phonon-assisted tunneling, we first de- 
fine these Dirac brackets as operators: 

c L = l°> n >> c Un = (5) 

then rewrite the mapping Hamiltonian in the operator 
representation: 

H cB = e Vkn{t)c\ kn C nkn +Yy e dn{t)c\n C dn 

-AV^+i(4„ +1 c dn + H.c.)] 

+ E^«( c W* c * 1 + ILc -)> ( 6 ) 

77, k : n 

with e n kn = ^k{t) + nhio ph and e dn = e d (t) + nfuv ph 



(n > 0). Note that these pseudo-Fermi operators in the 
effective Hamiltonian are bilinear. Finally, the explicit 
anticommunicators of these operators can be easily de- 
termined: 

{ c dn' 4ml = $ntni {Oqkn, c l'fc' m } = S nv ' 5kk' S nm - (7) 



B. Mathematical Formulation 

In the subsection, we derive the general formula of 
the time-dependent transmission T to t(t) and the time- 
average transmission Ttot and current / through the 
QD coupled to dispersionless Einstein phonons with fre- 
quency ujph by using the NGF technique^ The time- 
dependent current Jhn{i) from the nth channel of the left 
lead to the QD can be calculated from the time evolution 
of the occupation number operator of the electrons in this 
channel, N Ln = J2k c Lkn c Lkn' and one rea dily finds 

J Ln (t) = ~e(N Ln ) = - ([H eS (t),N Ln ]) 
n 

= T^ V ^,Lkn^t)-V£ n (Gt kn Jt,t))]- (8) 
k 

Here we define these hybrid retarded (advance) 
and lesser (greater) GFs G r n { ^ k < {>) {t,t') = 
((c d n(t)\cl k Jt<W^<(>) as follwos: G r n %Jt,t>) = 
±i6(±t T t')({c in (t), c\ km (f)}>, G< vkm (t,t>) = 
i(ikm(t')c d n(t)) ^d G> tVkm (t,t') = -i{c dn {t)cl km (t')); 
and similar definitions for ^im.n = 
((c v k m (t)\c dn (t'))) r ^<>'>. With the help of the 
Langreth analytic continuation rules^ these hy- 
bird GFs can be related to the GFs of the QD 

G± )l<(>) (*,*') = ((c dm (i)|cL(i')» r(a) ' <(>) , where 
G±\t,t') = ±iO{±tTt'){{c dm {t),c dn {t')}) 1 G< n (t,t') = 

i(4n(* / )c d ™(*)), and G> n (t,t') = -i{c dm {t)c\ n {t')). 
Then the time-dependent current JLn(t) becomes 



e r* 

J Ln (t) = t\Vlu\ 2 J2 / M G »»(Mi)0£^(*i.*) + ^„(Mi)0^ 
n k J ~°° 

(9) 



in which gJ~ kn (t,t') are the free-electron GFs in the 
nth pseudo-channel of the rj lead without coupling to the 
central region. 13 Employing the property G r —G a = G > ~ 
G < and Noticing the definitions of these GFs, Eq. @ can 
be simplified as 

J Ln (t) = |l^n| 2 E / dt l{ G nn(t,h)9Z kn {tl,t) 
n 1. J -00 



-9$ kn (t, h)G< n (h,t)} . (10) 
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In the wide band limit, where the hopping matrix element 
Vnn — Kj is independent on energy, one has 

k J 

ElVI 2 ^»(*^') = -<r, / f^[i-/>)]F(M'), 

With F(t,t') = e -Mt-t')-iJ t Urv n (r)^ r?? 

2tt J2k l^?! 2 *^ — e »?fe) being the generalized linewidth 
function, and /" = {1 + e ( u +"H».-M,)/teTj-i being the 
Fermi distribution function of the nth pseudo-channel in 
the r\ lead. 

Therefore to compute the time-dependent current, we 
have to evaluate the Keldysh GFs G„n (t,~t')- In the 
following, we first calculate the retarded GF G r mn (t, t') in 
the equation of motion approach. Introducing the gauge 
transformation's 

G r mn (t,t>) = e->ti d ™^G r mn (t,t>), (12a) 
G r nkm Jt,t>) = e- l ti d ™^G r nkm Jt,t>), (12b) 

one easily finds 

~ e<2 - mfkuph^j G r mn (t, t') = S mn 5(t - t') 

-\Vrn + lG r (m+l)n (t,t') - \<SmG\ m _ 1)n (t,t') 
+ Y / V vm G r vkmtn (t,tiyK d ^\ (13) 

rj . k 



and 



Q \ _ _ 



with A,,(i) = Vd(t) — v n (t). Then we transform the time 
coordinates to energy coordinates following the usual pre- 
scription, 

Y(u, uf) = / dtj dt' Y(t, t'y^-^'t' , (15) 
and obtain 

[u-ed - mhu ph - T? m {uS)] G r mn {u), u)') = 5 mn S(uj - u') 
-\y/mG r {m _ 1)n {u,u') - \Vm + lG[ m+1)n (w, a/), (16) 

where the retarded self-energy £^(0;) due to the coupling 
to the leads through the same gauge transformation as 
for the GFs is defined as 

= Ei^wi 2 ^*m( w ' w, )*(w-«') 

77, k 



E 

77, fc 



iK, m | 2 



(17) 



In the wide band approximation, we have = — | (r^ + 
r R ) = -|r, which leads to — > -i§5(f - *')■ 
Here to avoid unnecessary complications, we ignore the 
EPI induced band narrowing, which will lead to some 
trivial quantitative but no significant qualitative changes 
for tunneling current^ Therefore, it is helpful to rewrite 
the resulting Eq. (|16fl in a compact matrix form: 



[{w + iT)I -B]G r (Lu) = I 1 (18) 

in which B is a TV x N (N = 00 denotes the total phonon 
number) symmetrical tridiagonal matrix: B nn — e^+in— 
l)hujph, Bn(n-i) = -\y/n - 1, and S„ ( „ +1) = -A\/n, I 
is the N dimensional unit matrix. A similar result is ob- 
tained for the advanced GF G a (ui), and simply we have 
G a = [G r ]t in the Fourier space. After performing the 
reverse Fourier transformation of G r (uj), then substitut- 
ing the resulting G r (t, t') into Eq. I|12al) . one obtains 



G r (t,t') 



du> 



1 



2tt (u + iT)I - B 



-ioj(t— t')— i J,/ drvd{r) 



(19) 

Now we proceed to calculate the correlation GFs 
G< >s >(t,t'). Using the formal Keldysh GF technique, 
they are related to the retarded and advanced GFs as 

G< > \t,t') = Jdh y^ 2 G r (t,i 1 )S<(»(t 1 ,i 2 )G a (t 2 ,0, 

(20) 

with the "scattering in (out)" Keldysh self-energy 
in the matrix form with respect to the phonon num- 
ber. Because the strongly correlated Hamiltonian is 
transformed exactly to an one-body problem these 
Keldysh self-energies are produced only by tunneling cou- 
pling to the leads and they are easily derived, in the wide 
band limit, as 



n 

= E E £>)*(*-*%™ 

n 

E^Il-^Ml^t-fMrnn 



y> 

ran 



(t,0 



(21a) 



(21b) 



With Eqs. US) and (J23J, the Keldysh GF G^> (t, t') can 
be determined 

G <{>) (t t') 



Ei du 

x^J^^^tJ^^fOr, (22) 
in terms of the the auxiliary function AL mn (u,t): 



—iu>(t — t')—i J.j drv v {r) 



dti 



du' 
2T 



J ( W '-w)(t-t 1 )-i/ t t i dr A„(t) 



(23) 
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Finally we substitute Eqs.JTTJ, (J2H and (J23) into 
Eq. ljTU)) to calculate the current Jl„(£). As described 
in the above subsection, the current in the nth pseudo- 
channel of the left lead can be divided into two parts: 
Jj~n(t) denotes the component outgoing from the Ln 

channel to all pseudo-channels in both leads while Jj^ (t) 
means the part incoming from all pseudo-channels in 
both leads into this counted channel. It is obvious that 
the two terms bracketed in a pair of parentheses in 
Eq. IjlUfl are exactly corresponding to the two different 
tunneling processes, and we will write them down sep- 
arately in the following. After some algebra, we obtain 



= l-r^E /^/£(^)S> i ^)|A; n/ ( W ,t)| 2 ,(24a) 

n ,l J 

(24b) 



'/■' 



Another important point we must emphasize is that an 
appropriate weight factor should be added for any com- 
ponents: product by the probability P belonging to the 
pseudo-channel electrons inject from. Therefore, the to- 
tal time-dependent current flowing through the left lead 
is as a sum over all pseudo-channels in the left leads: 



m = I^l J ^ek«^)i 



x {P„/£Hiyi - 4(c)] - - /£(«)]} . (25) 

Taking into account the fact that A^(u>, t) is a symmetric 
matrix, the total current J{t) is simplified as 



J(t) = lr L r R J d^Y.WRnl^M 



x {P„/£(u,)[l - f l R (u)] - Pif R (u>)[l ~ f2(u>)]} . (26a) 

We can also compute only the elastic contribution to the 
total current by imposing the constraint of elastic tun- 
neling n = I: 

Jel(t) = j-TlTr I d^\ A Rnn(u,t)\ 2 P n {f2(u)~fR(u)}- 
J n 

(26b) 

Similarly, we can define the time-dependent total trans- 
mission Ttot(w,t) and the elastic transmission T e i(u>,t) 
as: 



T tot (w,t) = T L T R ^2P n \A Rnl (u,t)\ 2 , (27a) 

n,l 

T d (u,t) = T L T R ^2P n \A Rnn (uj,t)\ 2 . (27b) 

n 

In order to check our model and theory, we inspect the 
derived current and transmission probability formulae in 
two special cases: 1) without the MW irradiation field; 



2) absence of the EPI. First if all time-varying parts in 
energy are zero v LjRA = 0, one has A r vnl (u),t) = G r nl (u). 

Then t n i(uj) = T lY R \G r rd (w)\ 2 describes the transmission 
probability from nth channel to the Zth channel accord- 
ing to the Lee-Fisher formulation^ The total transmis- 
sion becomes T tot (uj) = ^2 n iP{n)t n i{uj), and the time- 
independent total current is 



J 



n , I 



x {p n /£H[i - /AM] - Pifk(o>)[i - f2(<»)]} ■ (28) 

Both of them are in consistent with the previous results 
based on the Landauer-Biittiker formulation. Further, 
if there is no EPI, A = 0, the matrix B reduces to a 
diagonal matrix with the nonzero element B nn = Ed + 
(n — l)hw p h. Easily, we obtain 



J=Jc 



-T T T 



h L R J ^ | lu — td — nhjjph + iT | 2 
xP n [/rM-/jj(u,)]. (29) 



1 



Recalling J2 n P n = 1, simple algebra calculation gives 

1 



J = T T L T R 
h 



dio- 



e d + iT 



-/«(")]■ (30) 



This is the famous time-independent current formula for 
the non-interacting resonant-level model. 

We return to the photon-phonon assisted current cal- 
culation. Because the easily measured current is the dc 
component, we must compute the time-average current 
I = (J(t)). Define the time-average of a time-dependent 
object Y(t) as 



(Y(t)) 



- Lim — 

t p * OO £ , 



tp/2 



dtY(t). 



(31) 



P J-t p /2 



Of course, if Y{t) is a periodic function of time, it is 
sufficient to average over the period. For the harmonic 
time modulation interested here, the integration over t\ 
in the auxiliary function A^ mn (cj,t) is easily to carried 
out: 



exp 



A 

x E Jl(-£)e im G r mn (u,-m), (32) 

/— — oo 

where Ji(x) is the Ith order Bessel function of the first 
kind and A,, = Vd — v v is the difference of the intensities 
of the MW field on the QD and the n lead. Therefore, 
the time averaging of the current and the transmission 
TtotM = (Ttot(uj,t)) will be 

/ OO A 
^E E j?(-^)\G r nm ("-im 2 



n,m [— — oo 
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x {P n fl{u)[l - f%{u>)\ - P TO /£M[1 - /£(*)]} , (33) 
and 



T tot ( w ) = r i r fl ^p„ £ J?(-jf)|<5; m (w-/n)| 2 , 

n,m i= — oo 

(34) 

and two similar expressions for I c \ — (J c \(t)) andT i(co>) = 
(T c i(w,t)) if setting n = m. 



III. NUMERICAL RESULTS AND 
DISCUSSIONS 

In this section, we start to numerically study the 
properties of the time-average transmission and current 
through a single-molecular QD under MW irradiation 
fields, on the basis of Eqs. I|33[) and i|34|) . For simplic- 
ity, in our calculation we assume the tunneling coupling 
between the molecular QD and the two leads to be sym- 
metric Tl = Tr = r, and the applied MW fields to be 
also symmetric vl = vr = v and then A^ = A# = A. 
We set the phonon energy as the unit of energy through- 
out the rest of the paper and choose the Fermi levels of 
the two leads as the reference of energy /iL = fiR = at 
equilibrium condition. 

A considerably wide range for the values of the phonon 
energy has been estimated for various mesoscopic sys- 
tems in some experimental papers, from 10 fieV to 
lOmeV, in addition to a very weak coupling to the leads 
r ~ 1 fxeY In the literature, some different val- 
ues for the phonon energy u p h, EPI constant A, and 
the tunneling constant T have been used for theoretical 
analysis i2*2i2iiliiiiiiLi& In the present paper, we will par- 
ticularly choose these parameters as typical values for the 
purpose of providing a satisfactory interpretation for the 
obtainable experimental results about low-temperature 
tunneling in the single-molecular QD. We choose A = 0.5, 
r = 0.04, and aj p h — 1 as the unit of energy. In the fol- 
lowing calculations, the temperature is assumed to be 
T = 0.04. 

Worth pointing out that the present approach is an 
exact method to deal with the electron-phonon coupled 
system with arbitrary coupling strength if the maximum 
number of phonons N p h involved in calculation is counted 
up to infinite. Unfortunately it is impossible in real cal- 
culation and we have to truncate this number to a finite 
value as long as the computation convergence is reached 
with any desired accuracy. The appropriate value of N p h 
depends on the energy of the Einstein phonon mode, the 
EPI constant, and the temperature of the system under 
investigation. Concerned with these parameters in the 
present paper, we choose N p h = 5 to obtain results with 
higher than 1% accuracy. 



A. Time-average transmission 

First we examine the time-average transmission. In 
Fig. 2 we plot the total and clastic time-average trans- 
missions as a function of the incoming electron energy 
u) for the case of QD bare level e<j = under different 
MW irradiation fields: the frequency £1 of the MW field 
is (b) lower than, (c) equal to, and (d) higher than the 
frequency of the dispersionless phonon w p h for the case 
of a weak MW intensity A/fi < 1 (A = 0.2). For com- 
parison, we also plot the results without irradiation in 
Fig. 2(a). The pair of numbers (n, m) located near the 
peaks means numbers of the phonon and photon elec- 
trons exchange with the heat bath and the MW field, 
respectively, when electrons tunnel through the central 
region. The positive (negative) of the number denotes 
emission (absorption) of quanta. In absence of irradia- 
tion fields, three important features can be observed with 
respect to the EPI: 1) The central peak of the transmis- 
sion shifts from the position ui = to e c i — A 2 /co p h with 
a slightly suppressed amplitude (smaller than 1) due to 
the polaron effect; 2) The positions of the side peaks cor- 
responding to emission of n phonons are approximately 
given by u> = e^ — A 2 /ui p h + nHu) p h, meaning that the new 
pseudo-channel is opened and participates contribution 
to tunneling; 3) There are no peaks in the left side of the 
central peak because before tunneling there is no avail- 
able phonon for absorption at low temperature. These 
observations are just in agreement with the previous the- 
oretical results£i2iiLi& 
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FIG. 2: Time-average transmission probability as a function 
of the incident electron energy under different irradiations: 
(a) No irradiation; (b) The frequency of the MW field Q = 0.5; 
(c) 0=1; (d) Q, = 1.5. The intensity of the irradiation is 
A = 0.4. The solid line is the total transmission T to t, and 
the dashed line is the elastic part T e i. The pair of numbers 
(n, m) near each peak denotes the numbers for emission (pos- 
itive number) or absorption (negative number) of phonon and 
photon, respectively, involved in the tunneling process. The 
other parameters used in calculation are: V = 0.04, td = 0, 
Lo ph — 1.0, and A = 0.5. The temperature is set to T = 0.04. 
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In the presence of irradiation MW field, it is clear to 
observe that more rich peaks have been observed in the 
transmission when the involved quanta pair (n, m) sat- 
isfies uj — €d — A 2 /ujph + nftiOph + mhQ, for the given 
frequency of MW field. In contrast to the phonon, pure 
absorption peaks of photon can be seen in negative in- 
cident energy region due to the fact that photon is in- 
depdendent of temperature. For the case of lower MW 
frequency f2 = 0.5 < uj p h [Fig. 2(b)], there appears a new 
peak in the middle of the zero-phonon-peak (0, 0) and the 
first satellite-phonon-peak (1,0), which is attributed to 
the following two processes: (0, 1) pure emission of one- 
photon and (1,-1) emission of phonon assisted by ab- 
sorption of one-photon. And these phonon peaks, such as 
(1,0) and (2,0), could be also accompanied by emission 
of photon and photon-assisted multi-phonon processes, 
leading to obvious enhancement of the contribution of 
the elastic channel to transmission, as shown in Fig. 2(b) 
and (c) for the cases of irradiation fields with low fre- 
quencies. If the frequency of the MW field is higher than 
the frequency of phonon, no photon-assisted processes 
take place between the main phonon peak (0, 0) and the 
first phonon peak (1,0) because there is insufficient en- 
ergy to excite the high-frequency photon. As a result, the 
contribution of the elastic channel keeps low. More in- 
terestingly, we observe a weak emission of one-phonon at 
the negative energy region with the help of absorption of 
one-photon even under this small intensity of irradiation. 




FIG. 3: Time-average transmission probability as a function 
of the incident electron energy under irradiations with fre- 
quencies lower [(a, b): £1 = 0.5] and higher [(c, d): fl = 1.5] 
than the frequency of the Einstein phonon. The intensities of 
the irradiations are A = 0.8 (a, c) and 1.6 (b, d). The other 
parameters are the same as in Fig. 2. 

Increasing intensity of the MW field will enhance the 
photon-assited processes. In Fig. 3 we plot the time- 
average transmissions under the MW fields of frequencies 
fl = 0.5 and 1.5 with stronger intensities A = 0.8 and 



1.6. Obviously, we observe more peaks with more high 
amplitudes involving emission and absorption of photons, 
as well as suppression of the pure phonon-assisted pro- 
cesses. For the strongest intensity A = 1.6 calculated 
here, the high-frequency photon (f2 = 1.5) even excites 
an emission peak involving two-phonon process between 
the main-phonon and the one-phonon peaks, as depicted 
in Fig. 3(d). 



B. 



Time-average current and nonlinear differential 
conductance 




FIG. 4: The time-average total current I (solid line), elas- 
tic current I c \ (dashed line) [(a-d)], and the corresponding 
differential conductance dl/dV [(e-h)] as a function of the 
applied bias voltage for a single-molecular QD with bare level 
ed = 1.0 under different irradiations: (a,e) No irradiation; 
(b,f) The frequency of the MW field Q = 0.5; (c,g) 0=1; 
(d,h) Q = 1.5. The intensity of the irradiation is A = 0.4. 
The other parameters are the same as in Fig. 2. 

According to the current formula Eq. (|33|l , nonlinear 
differential conductance, defining as the derivative of the 
time-average current with respect to the bias voltage 
dl/dV (dI e \/dV), is believed to be a possible tool in 
experiments to detect the photon-phonon-assisted multi- 
peaks in time-average transmission due to its proportion- 
ality to T tot ( e n(eV) at low temperature. So we illustrate 
in Fig. 4 the calculated time-average current [(a)-(d)] and 
the differential conductance [(e)-(h)] for the case of bare 
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level ed = 1.0 without ac field (a,e) and under ac irradi- 
ations with different frequencies SI = 0.5 (b,f), 1.0 (c,g), 
and 1.5 (d,h) at a weak intensity A = 0.4. For simplicity, 
we assume the bias voltage is added symmetrically, i.e., 
+V/2 on the left lead and -V/2 on the right lead. Ob- 
viously, the time-average current shows steps indicating 
that new satellite-peak in transmission denoted by (n, m) 
participate in conducting electrons. And correspond- 
ingly, there appears a peak in the dl/dV-V curve located 
at the bias voltage eV/2 = e^ — A 2 /uj p h + nhu> p h + mft£L. 
For the low frequency irradiation field SI = 0.5, we ob- 
serve a peak below the main-phonon peak (0, 0) resulted 
from the pure one-photon absorption process (0,-1), 
which is nearly an elastic phonon process and has lit- 
tle inelastic contribution to the current. Moreover there 
is another peak between the main-phonon and the one- 
phonon peaks, which can be attributed to the combina- 
tion contributions of a pure one-photon emission process 
(0, 1) and a one-photon-absorption-assisted one-phonon 
emission process (1,-1). It possess of course inelastic 
contribution with a small portion. 

As pointed above, in order to observe more photon- 
phonon-assisted processes, a possible method is to in- 
crease the intensity of ac field. In Fig. 5, we show the 
bias voltage-dependent time-average current and the dif- 
ferential conductance under several different microwave 
intensities A = 0.4, 0.8, and 1.6 for SI = 0.5 (a,b) and 
1.5 (c,d). We find the following features in this figure: 
1) Increasing irradiation intensity decreases the time- 
average current at higher bias voltage, but enhances the 
current at lower bias voltage; 2) The peak due to the 
pure phonon process (0,0) is suppressed while the peak 
concerned with the photon-sbsorption processes (0,-1) 
[Fig. 5(b)] or (1,-1) [(d)] grows up gradually with ris- 
ing of ac intensity; 3) These peaks due to photon-phonon 
combination processes have complicated relation with the 
ac intensity in the case of low-frequency irradiation field, 
which can be attributed to summation of the Bessel func- 
tion with respect to the ac intensity in the time-average 
current formula Eq. J33J). 

Now we turn to study the variation of time-average 
current with the irradiation frequency at a given bias. 
Fig. 6 depicts this MW spectroscopy for the QD with bare 
level ed = 1.0 at a low bias voltage eV = 0.2. In this case, 
the low-frequency irradiation field induces a remarkable 
enhancement of the time-average current, in comparison 
with the dc current Iq without irradiation, when elec- 
trons absorb energy of photon to overcome the energy 
gap between the left lead and the QD. The enhanced dc 
current is called pumped current. 2 — In Fig. 6, we observe 
two obvious peaks with one order of magnitude higher 
than the dc current Iq , which are due to the one-photon- 
and two-photon-absorption induced pumping. But the 
positions of the two peaks marked by triangle symbols ex- 
perience blue-shift SI = (eV/2+X 2 /iVpi l —ed)/m (SI ~ 0.71 
for m — 1, and SI ~ 0.35 for m — 2) due to the po- 
laron effect. More interestingly, we find two more ev- 
ident pumped peaks (1,-1) and (2,-1) with smaller 
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FIG. 5: The time-average current I [(a,c)], and the corre- 
sponding differential conductance dl/dV [(b,d)] vs bias volt- 
age for the case of = 1.0 under several different irradiation 
intensities A = 0.4 (dotted line), 0.8 (dashed line), and 1.6 
(solid line). (a,b) are ploted for SI = 0.5, (c,d) for Q, = 1.5. 



amplitude at high-frequency region, where the irradi- 
ation field can not only produce pumped current but 
also provide energy for electrons to emit phonons when 
electrons are pumped to tunnel through the molecule. 
These phonon-assisted pumped peaks locate at il — 
(eV/2 + X 2 /oj p h — ed~nujph)/'m. And they have of course 
higher amplitudes in the presence of stronger irradiation. 
This addresses that the MW spectroscopy of the pumped 
current is another possible method to observe the photon- 
phonon- assisted resonant tunneling. 
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FIG. 6: The renormalized time-average total current / vs 
the frequency of irradiation ac field for ed = 1.0 at the bias 
eV = 0.2 under several different irradiation intensities A = 
0.4 (solid line), 0.8 (dashed line), and 1.6 (dotted line). Io is 
the dc current of this system without irradiation. 
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IV. CONCLUSIONS 

In this work we study the low-temperature time- 
dependent resonant tunneling through a single- molecular 
QD subjected to a dispersionless local phonon mode with 
a strong coupling to electrons. Due to the weak elastic 
parameters in molecular materials, the energy of boson 
field is matched with the MW irradiation field. So when 
electrons enter into the tunneling region, electrons can 
easily excite the phonon field with and without assis- 
tance of the ac field. Therefore, it is predicted that the 
combination influence of the phonon and photon fields 
will essentially change transport properties in the sigle- 
molecular QD. 

Owing to the strong EPI, the traditional perturbation 
theory is invalid for this problem even though it has been 
intensively applied for interpretation the classic trans- 
port phenomena in bulk semiconductor. Following the 
pioneer work of Bonca and Trugman, we first transform 
this many-body problem exactly into a one-body scat- 
tering problem, by projecting the original Hamiltonian 
in the representation of electron-phonon coupled Fock 
space, i.e., the direct-product states of electron states 
and phonon number states. Based on this noninteract- 
ing Hamiltonian, we use the NGF approach to derive the 
time-dependent and the time-average transmission and 



current in the wide band approximation. For the sake of 
simplicity, we neglect the EPI-induced band narrowing in 
the present calculation. Moreover we adopt the adiabatic 
approximation for the irradiation MW field. 

Our numerical results show that time-average trans- 
mission vs incident electron energy u> displays additional 
peaks, besides original phonon absorption peaks, due 
to photon emission or absorption mediated processes as 
long as the irradiation frequency matches the condition 
u> = ed — A 2 /u)ph + nhujph + mhfl. We also point out 
that these new features will be observed in experiments 
by measuring the nonlinear differential conductance and 
the MW spectroscopy of pumped current at proper irra- 
diation intensity. 
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